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Abstract 

We advocate charmed-hadron inclusive hadroproduction as a laboratory to probe intrinsic charm 
(IC) inside the colliding hadrons. Working at next-to-leading order in the general-mass variable- 
flavor-number scheme endowed with non-perturbative fragmentation functions recently extracted 
from a global fit to e~^e~ annihilation data from KEKB, CESR, and LEPl, we first assess the 
sensitivity of Tevatron data of D^, D^, and D*^ inclusive production to the IC parameterizations 
provided by Pumplin et al. We then argue that similar data from pp collisions at RHIC would have 
the potential to discriminate between different IC models provided they reach out to sufficiently 
large values of transverse momentum. 
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I. INTRODUCTION 



Recently, the inclusive production of charmed hadrons (Xc) at hadron colliders has been 
the subject of extensive experimental and theoretical studies. The CDF Collaboration mea- 
sured the differential cross section da/dpx for the inclusive production of D^, D^, D*^, 
and mesons (and their antiparticles) in pp collisions at the Fermilab Tevatron (run II), 
with center-of-mass energy ^/s = 1.96 TeV, as a function of the transverse momentum (p^) 
in the central rapidity (y) region [l]. The PHENIX Collaboration measured non-photonic 
electron production through charm and bottom decays in pp, dAn, and AuAu collisions at 
the BNL Relativistic Heavy Ion Collider (RHIC) with ^/s = 200 GeV as a function of pt Q] ■ 
The STAR Collaboration at RHIC presented mid-rapidity open charm spectra from direct 
reconstruction of /if* K^n^ decays in dAu collisions and indirect electron-positron 
measurements via charm semileptonic decays in pp and (iAu collisions at ^/s = 200 GeV 
0. Recently, they also reported result, on non-photonic electron ptodnctton in pp, .An, 
and AuAu collisions at y/s = 200 GeV [3|. Unfortunately, these RHIC data only cover a 
very limited small-pr range, where theoretical predictions based on perturbative QCD are 
difficult. 

On the theoretical side, the cross sections for the inclusive production of X^. mesons can 
be obtained as convolutions of universal parton distribution functions (PDFs) and universal 
fragmentation functions (FFs) with perturbatively calculable hard-scattering cross sections. 
The non-perturbative input in form of PDFs and FFs must be known from fits to other 
processes. The universality of the PDFs and FFs guarantees unique predictions for the 
cross section of the inclusive production of heavy-fiavored hadrons. The results of such 



calculations for production at t 



le Tevatron have been presented recently by us in Ref. 



and compared to the CDF data Ij. For all four meson species, D^, D^, D*^, and Df, 
we found good agreement with the data in the sense that the experimental and theoretical 
errors overlap. For the D^, D*~^, and mesons, many of the central data points fall into 
the theoretical error band due to scale variation. Only the central data points for the 
mesons lie somewhat above it. For the D^, D^, and D*~^ mesons, the central data points 
tend to overshoot the central theoretical prediction by a factor of about 1.5 at the lower end 
of the considered px range. With the exception of the case, the experimental results are 
gathered on the upper site of the theoretical error band, corresponding to a small value of 
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the renormalization scale and large values of the factorization scales fip and /i^, related 
to the initial and final states, respectively. 

In our analysis Sj, we employed the general-mass variable-flavor-number scheme (GM- 
VFNS), which combines the zero-mass variable-fiavor-number scheme (ZM-VFNS) and the 
fixed-fiavor-number scheme (FFNS). The GM-VFNS is close to the ZM-VFNS, but retains 
iTi^/pr power terms in the hard-scattering cross sections. Here, m stands for the mass of the 
charm quark. The calculational details for hadron-hadron collisions were elaborated in our 
previous works j^. The characteristic feature of the GM-VFNS is that the charm quark is 
also treated as an incoming parton originating from the (anti)proton, leading to additional 
contributions besides those from the light u, d, and s quarks and the gluon. This is quite 
analogous to the ZM-VFNS, but with the important difference that the power terms are 
retained in the hard-scattering cross sections. This allows us to apply the GM-VFNS also 
in the region of intermediate pt values, pt ^ rn. In fact, in Ref. [sl], we had pt ^ rn. 

In Ref. jsl, we included n/ = 4 active quark flavors and took the charm-quark mass to 
be m = 1.5 GeV. The strong-coupling constant was evaluated in next-to-leading 

order (NLO) with A^^^ = 328 MeV, corresponding to a^^\mz) = 0.1181. We employed 
proton PDF set CTEQ6.1M from the Coordinated Theoretical-Experimental Project on 
QCD (CTEQ) Collaboration [?] and the FFs from Ref. |8i]. The default choice for the 
renormalization and the initial- and final-state factorization scales was fiR = fip = f^'p = "^t, 
where rriT = \/pt ~^ is the transverse mass. 

Compared to our previous work Sj, we have now much more reliable FFs for the tran- 
sitions u, d, s,c,g Xc at our disposal. The FFs used in Ref. ^ were determined by 
fitting_^the fractional-energy spectra of the hadrons measured by the OPAL Collabora- 



tion 
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m e+e~ annihilation on the Z-boson resonance at the CERN LEPl collider. 
These data had rather large experimental errors and the disadvantage of being at a rather 
large scale /i^ = Mz, far away from the scales of Xc production in pp collisions at the Teva- 
tron, which typically correspond to pt values below 25 GeV. In the meantime, new data on 
charmed-meson production with much higher accuracy were presented by the CLEO Collab- 



oraton 



11| at CESR and by the Belle Collaboration 12|] at the KEK collider for B physics 



(KEKB). These data offered us the opportunity to determine the non-perturbative initial 
condition of the FFs much more accurately. The data from CLEO 11] and Belle [l^ are 
located much closer to the threshold y/s = 2m of the transition c — >• than those from 
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OPAL |9|, 110|. The scale of the CLEO and Belle [12| data is set by the cm. energy 
-v/i = 10.5 GeV. Recently, Kneesch and three of us extracted from these data FFs for D^, 

n 

D^, and D*^ mesons; the details may be found in Ref. [13l |. 

Another important ingredient for the theoretical description of inclusive Xc production 
are the parton distribution functions (PDFs). Let fa{x,fiF) denote the PDF of parton a 
inside the proton at momentum fraction x and factorization scale fip. At short distances, 
corresponding to large values of fip, the scale dependence is determined by the QCD evo- 
lution equations with perturbatively calculable evolution kernels. So, the PDFs are fully 
determined by their functional forms in x specified at a fixed scale fip = f^o provided fiQ 
is large enough to be in the region where perturbative QCD is supposed to the valid. In 
most applications, /iq is chosen to be of order 1-2 GeV, which is at the borderline between 
the short-distance (perturbative) and long-distance (non-perturbative) regions. The PDFs 
of the gluon and the light quarks (a = g,u,d,s) are non-perturbative. They are obtained 
phenomenologically through global QCD analyses, in which the theoretical predictions are 
compared with a wide range of experimental data on hard processes 

m 



The cross sections of inclusive charmed-meson production in the GM-VFNS depend heav- 
ily on the PDF of the charm quark and less on those of the light quarks. In the following, we 
use the short-hand notation c{x, jj^p) = fc{x, /j^f)- In the global analyses of Refs. [3, Q, Q 
and many others jl^, the charm quark is considered a parton. It is characterized by the PDF 
c{x,fiF) that is defined for hf > rn. In common global QCD analyses at NLO in the MS 
scheme, the charm quark is considered as a heavy quark, and the ansatz c(x, /iq) = with 
fiQ = m is assumed as the initial condition for calculating c{x,fiF) at higher factorization 
scales /ii? > /iQ. This is usually referred to as the radiatively-generated-charm approach. In 
the global analyses cited above, this ansatz implies that the charm parton does not have any 
independent degrees of freedom in the parton parameter space, i. e., c(x, ^f) is perturbatively 
determined by the gluon and light-quark parton parameters. 

However, a purely perturbative treatment of the heavy quarks might not be adequate. 
This is even more true for the charm quark with mass m ^ 1.5 GeV, which is not much larger 
than typical hadronic scales of a few hundred MeV. In fact, there exists the possibility, that 
c{x, /io) 7^ 0. Actually, for many years, non-perturbative models exist th at g ive non-zero 



predictions for c(a;, /io) at some initial factorization scales fiQ = 0{m) [l7|, [l8|. The model 



of Ref. ITj was also analyzed in the framework of the FFNS 19|. In this connection, the 
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notion intrinsic charm (IC) has become customary. The IC models of Refs. |l7l . Il8[ | and a 
third one, which will be specified in a later section, were put to a stringent test by extending 
the recent CTEQ6.5 global analysis 1^ so that the charm sector has its own independent 
degrees of freedom at the initial factorization scale = m {2^. As a natural extension of 
the CTEQ6.5 analysis [14,], Pumplin et al. 20| determined the range of magnitude of IC 
that is consistent with an up-to-date global analysis of hard-scattering data with various 
assumptions on the shape of the x distribution of IC at a low factorization scale. 



The authors of Ref. 



20[ | found the IC of the light-cone models 



17 



to be compatible 



with the global data sample for magnitudes ranging from zero up to three times of what had 
been estimated in more model-dependent investigations. In these models, there can be a 
large enhancement of c{x,^p) at x > 0.1, relative to previous analyses that have no IC. The 
enhancement persists to rather large scales ^p-, up to 100 GeV. Therefore, these nonzero IC 
contributions can have an important effect on charm-initiated processes at HERA, RHIC, the 
Tevatron, and the LHC. Of course, at hadron colliders, the charm production cross section 
would be sensitive to c{x,iif) at a; > 0.1 only at sufficiently large values of xt = 2pt/\/s- 
For example, the cross section data at the Tevatron ij only cover the range pt < 25 GeV, 
so that Xt < 0.025, which is too small to be sensitive to IC contributions at x > 0.1. Only 
at RHIC, where ^/s is a factor of 10 smaller, we may expect sufficient sensitivity. At given 
values of ^/s and pt, the sensitivity to c{x, fip) at large values of x could be further enhanced 
by performing measurements at large values of \y\, i.e. in the extreme forward or backward 



regions 



2l|. 



The obvious data that could provide us with information on the IC contribution, would 
be those on the charm structure function F2 measured at HERA. But these data are already 
used in the global analysis [20!]. Unfortunately, they do not have a significant effect because 
of the rather large experimental errors and because the data are mostly located at small 
values of x. A comparison of such data with results from earlier CTEQ global analyses may 
be found in Ref. [22 1. 

After these general remarks, it is clear that it is interesting to study cross sections of 
inclusive charmed-meson production at RHIC {^/s = 200 GeV) and the Tevatron {^/s = 
1.96 TeV) more closely and to establish those kinematic regions which are sensitive to the 
IC component of the proton. For this purpose, we adopt the GM-VFNS approach outlined 



in Ref. 



We first compare the D^, D'^, and D*^ production data from Ref. [l| to 
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NLO predictions evaluated with the respective FF sets [13|] so as to test the latter (with or 
without IC). Then, we present results for the production cross section under Tevatron 
and RHIC experimental conditions assuming the six IC scenarios introduced in Ref. 20|. 

The outline of this paper is as follows. In Sec. [TTl we briefly review the FFs determined 
in Ref. [l^. In Sec. Illlt we describe the three IC models investigated in Ref. 20 1 as much 
as is needed to understand our final results. In Sec. HVl we present our NLO predictions 
for inclusive Xc production in pp and pp collisions. Section |V] contains a summary and an 
outlook. 



II. FRAGMENTATION FUNCTIONS FOR CHARMED MESONS 



For our calculation of the differential cross section d'^a/ {dpr dy) of p + p(j)) —^Xc + X, 
where Xc = D^, D'^, D*^ and X stands for the residual final state, a crucial ingredient is 
the non-perturbative FFs for the transitions a X^ where a = g,u,u,d,d, 3,3,6,15. For 



Xc = D*'^, such FFs were extracted at LO and NLO in the modified minima 
(MS) factorization scheme with Uf = 5 massless quark flavors several years ago 
distributions da/dx in scaled energy x = 2E / \f s of the cross sections of + e~ 



-subtraction 



23| from the 



measured by the ALEPH [2^] and OPAL [lO] Collaborations at LEPl. Two of us 



X 



25| 



extended the ^nalysis of Ref. (23[ to include Xc = D^,D~^,Df,A^ by fitting appropriate 
OPAL data [9]. Besides the total Xc yield, which receives contributions from Z —>■ cc 
and Z ^ bb decays as well as from light-quark and gluon fragmentation, the ALEPH and 
OPAL Collaborations separately specified the contributions due to tagged Z ^ bb events 
yielding Xf, hadrons, which then weakly decay to Xc hadrons. The contribution due to the 
fragmentation of primary charm quarks into Xc hadrons approximately corresponds to the 



difference o 
In Refs. 



these two measured distributions. 



3, 



25| . the starting point fiQ for the DGLAP evolution of the a Xc FFs in the 



factorization scale fi'p were taken to be /i = 2m with m = 1.5 GeV for a = g, u,u, d, d, 3,'s 
and /i = 2mfe with = 5 GeV for a = b,b. The FFs for a = g,u, u, d, d, 3, s were assumed 
to vanish at fi'p = fiQ and generated through DGLAP evolution to larger values of fi'p. For 



consistency with the MS prescription for the PDFs, these fits of the 



the choice /xq = m, nib in Ref. 



These determinations of Xr FFs 



Fs were repeated for 



23|, 



25| were all based 



solely on data from the Z-boson resonance. In these data, the effects of finite quark and 
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hadron masses were greatly suppressed and could safely be neglected. 

The most recent fits for Xc = D^, D^, D*^ reported in Ref. 13| include, beside the 



OPAL 



and ALEPH 



2J] data, much more precise data from CLEO ll[ and Belle 12|. 



They offered us the opportunity to further constrain the charmed-hadron FFs and test their 
scaling violations. However, the lower cm. energies of the CLEO and Belle [l^ data 
necessitates the incorporation of quark and hadron mass effects, which are then no longer 
negligible, into the formalism. The GM-VFNS, which is also the basis of the computation 
of the cross section of p + p Xg + X in Refs. js, Q|, provides the appropriate theoretical 



framework also for this (see Ref. [13| for details). We adopted the values of /zq and m from 
Ref. 3. 

In this framework, new FFs for and D*^ mesons were determined throu gh g lobal 



fits to all available e"'"e annihilation data, from Belle 121, CLEO llll, ALEPH 24fl, and 

nn 

OPAL [9|, 110|. In contrast to the situation at the Z-boson resonance, we had to take into 
account the effect of electromagnetic initial-state radiation on the Belle and CLEO data, 
which distorts the distribution in scaled momentum Xp = 2p/y^ of the cross section for 
continuum production in a non-negligible way [13]. We found that the global fits suffer from 
the fact that the Belle and CLEO data tend to drive the average of x of the c ^ X^ FFs to 
larger values, which leads to a somewhat worse description of the ALEPH and OPAL data. 
Since the 6 — > Xc FFs are only indirectly constrained by the Belle and CLEO data, namely 
via their contribution to the DGLAP evolution, their forms are only feebly affected by the 
inclusion of these data in the fits. In other words, the b ^ Xc FFs are essentially fixed by the 
ALEPH and OPAL data alone. It was found that hadron mass effects are more important 
than quark mass effects in the global fits. In fact, they are indispensable for fitting the lower 
tails of the Xp distributions from Belle and CLEO 

The z distributions of the c and b quark FFs at their starting scales were assumed to 



obey the Bowler ansatz 



26| 



(1) 



with three parameters N, a, and 7. Specifically, the fitting procedure was as follows. At the 
scale = m = 1.5 GeV, the c-quark FF was taken to be of the form of Eq. ([T]), while the 
FFs of the light quarks q {q = u, d, s) and the gluon were set to zero. Then these FFs were 
evolved to higher scales fi'p using the DGLAP equations with n/ = 4 active quark flavors 
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and =321 MeV. When the scale n'p reached the threshold value /i'^ = m;, = 5.0 GeV, 
the bottom flavor was activated and its FF was introduced in the Bowler form of Eq. ([1]). 
The evolution to higher scales n'p was then performed with nj = 5 and = 221 MeV. The 
values of the Bowler parameters Nc, ac, and 7c for c —>■ Xc and Nb, at,, and '-/b for b Xc 
thus obtained may be found in Tables 1, 2, and 3 of Ref. 13| for the D^, D^, and D*^ 
mesons, respectively, together with the achieved values of per d.o.f. (x^)- The values 
differ somewhat for the three D-meson species, but they are all acceptable. The smallest 
value was obtained for the meson. 



III. PDFS WITH INTRINSIC CHARM 



In a recent paper [20|, Pumplin et al. extended the up-to-date CTEQ6.5 global analysis 
14^ to include a charm sector with nonzero c(x, /io) at the initial factorization scale /iq = m. 
For this purpose, they considered three scenarios. The first two scenarios are based on the 
light-cone Fock-space picture of nucleon structure formulated in more detail by Brodsky . 
In this picture, IC is mainly present at large momentum fraction x, because states containing 
heavy quarks are suppressed according to their off-shell distance, which is proportional to 
[p^ + m?) / X. Thus, components with large m appear preferentially at large x. A wide variety 



of light-cone models predict similar shapes in x, as has been shown recently 28|]. Specifically, 
the light-cone models considered in Ref. [2^ include the original model of Brodsky, Hoyer, 
Peterson, and Sakai (BHPS) 171] and the so-called meson-cloud picture ISj, in which the IC 
arises from virtual low- mass meson-plus-baryon components of the proton, such as D°A+. 
The meson-cloud model also predicts a difference between c{x,fi'p) and c{x,fj,'p) and so 
provides an estimate of the possible charm/anticharm difference. 

Unfortunately, the light-cone formalism does not allow us to calculate the normalization 
of the uudcc components 28|, although estimates of the order of 1% were reported in the 



literature 



18 



29|. Therefore, in Ref. 20|, the magnitude of IC is determined by comparison 



with data incorporated in the global fit. 

The X dependence originating in the BHPS model 



IZl is 



c(a;, /xo) = c(a;, /xq) = v4a;^[6x(l -|- x) Inx + (1 — a;)(l + lOx -|- x^ 



(2) 



In Eq. ([2]), /io = m with m = 1.3 GeV and the normalization constant A is treated as a free 
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parameter in the fit. Its magnitude is controlled by the average c + c momentum fraction, 



(x) 



c+c 



dxa;[c(x,/io) + c(a;,/io)]. 



(3) 



The exact x dependence predicted by the meson-cloud model cannot be given by a simple 
formula. However, it was shown in Ref. 



28l | that the charm distributions in this model can 



be very well approximated by 



(4) 



where the normalization constants A and A are constrained by the quark number sum rule. 



dx [c(x,/io) - c(x,/io)] = 0, 



(5) 



which fixes A/ A. The overall magnitude of IC is again determined in the global analysis. 

As a third scenario, Pumplin et al. 20| studied a purely phenomenological scenario in 
which the shape of the charm distribution is sea-like, i.e. similar to those of the light- 
flavor sea quarks, except for an overall mass suppression. In particular, they assumed that 
c(a;, jUo) = c(x, /io) oc /iq) + c?(x, /iq) at ji^ = m. 

The c(a;, /io) and c(x,/io) functions of the BHPS, meson-cloud, and sea-like models dis- 



cussed above are used in Ref. 
as explained in Ref. 



20| as input for the general-mass perturbative QCD evolution 



141 ]. Then, the range of the IC magnitude is determined to be consis- 



tent with the global data fit. The quality of each global fit is measured by Xgiobai; which is 



shown in Fig. 1 of Ref. |20] as a function of {x)c+c for the three models. From this figure, 
one can see that in the lower range, < {x)c+c ^ 0.01, Xgiobai varies very little, i.e. the 
fit is very insensitive to {x)c+c in this interval. This means that the global analysis of the 
hard-scattering data gives no evidence either for or against IC up to {x)c+c ~ 0.01. For 
{x)c+c > 0.01, the curves for the three models in Fig. 1 of Ref. |20| rise steeply with {x)c+c- 
Xgiobai ~ 3450 represents a marginal fit in each model, beyond which the quality of the fit 
becomes unacceptable according to the procedure established in Refs. j7,[l^. This criterion 
is based on the fact that one or more of the individual experiments in the global fit is no 
longer fitted within the 90% confidence level. For comparison, the fit with no IC yields 
Xgiobai ~ 3330. The {x)c+c distributions of Xgiobai achieved in the BHPS, meson-cloud, and 
sea-like models reach the marginal-fit limit of 3450 at 2.0%, 1.8%, and 2.4%, respectively. 
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The authors of Ref. 



201 ] studied also the typical, more moderate {x)c+c values 0.57%, 0.96%, 



and 1.1% in the same order. The x dependences of c{x,^f) and c{x,^f) obtained in the 
three models are shown for the two values of {x)c+c mentioned above and for various val- 
ues of fip in Figs. 2-4 of Ref. 2o| and are compared there with the zero-IC result of the 
CTEQ6.5 fit. From these figures one notices that, in the two light-cone models, IC leads 
to an enhancement of c{x,^f) and c{x,^f) at x > 0.1 relative to the PDF analysis with 
zero IC, while the deviation is small for a; < 0.1. In the sea-like model, there is also a 
significant enhancement of c{x,^f) and c{x,^f) relative to the zero-IC PDF. In this case, 
the enhancement is spread more broadly in x, roughly over the region 0.01 < x < 0.50. 



le three models with two values 
, are available together with 



20 



16| as seven members, which 



The six PDF sets constructed with IC, corresponding to t 
of {x)c+c each, which are designated by CTEQ6.5C in Ref. 
the corresponding zero-IC PDF set via the LHAPDF standard 
we denote as CTEQ6.5Cn with n = 0, 1,2,3,4,5,6. Specifically, n = stands for zero IC, 
n = 1, 2 for the BHPS model with {x)c+c = 0.57%, 2.0%, n = 3, 4 for the meson-cloud model 
with {x)c+c = 0.96%, 1.8%, and n = 5, 6 for the sea-like model with {x)c+c = 1-1%, 2.4%. 



IV. COMPARISON WITH CDF DATA AND PREDICTIONS FOR RHIC 
A. Predictions for CDF 

In this section, we present our predictions for the differential cross section d^cr/(dpT dy) 
oi p + p ^ + X with Xr, = D^,D+, D*+ at NLO in the GM-VFNS. The calculation 
proceeds as outlined in Ref. [gj with the modifications explained in Sec. HIl 

The NLO cross section consists of three classes of contributions. 

1. Class (i) contains all the partonic subprocesses with a c, c — > X^ transition in the final 
state that have only light partons {g,q,q) in the initial state, the possible pairings 
being gg, gq, gq, and qq. 

2. Class (ii) contains all the partonic subprocesses with a c, c — transition in the final 
state that also have c or c quarks in the initial state, the possible pairings being gc, 
gc, qc, qc, qc, qc, and cc. 

3. Class (iii) contains all the partonic subprocesses with a. g,q,q ^ Xc transition in the 
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final state. 

The contributions of classes (ii) and (iii) are calculated in the ZM-VFNS using the hard- 
scattering cross sections usually used for the inclusive production of light mesons 30| . The 
light-quark fragmentation contributions are negligible. However, gluon fragmentation con- 
tributes significantly, as was first noticed in our previous work jsj. 



As was shown in Ref. [g], the m dependence of the class- (i) contribution is greatly screened 
in the full cross section because the latter is dominated by the contributions of classes (ii) 
and (iii), which are calculated with m = 0. In fact, the bulk of the contribution is due to 
class (ii), which provides a handle on the c(x, /i^) and c{x,fiF) and thus on IC. 

The CDF data [ll come as distributions da/dpT at ^/s = 1.96 TeV with y integrated 
over the range \y\ < 1. For each Xc meson, the particle and antiparticle contributions 
are averaged. We work in the GM-VFNS with n/ = 4 thus excluding Xc hadrons from 
Xft-hadron decays, which are vetoed in the CDF analysis 1|. For the comparison with the 



CDF data, we use set CTEQ6.5 of proton PDFs [ij] without IC (see Sec. IIIip . and set 



global-GM of Xc FFs ll3| (see Sec. We adopt the value = 328 MeV, which yields 



af'\Mz) = 0.1181, from Ref. [ij] and the values m = 1.5 GeV and = 5 GeV from 



MS 



Ref. 
Ref. 



There is a slight mismatch of the input parameters Uf, A 



(4) 



MS- 



141 ] uses Uf = 5, m = 1.3 GeV and nib = 4.5 GeV and Ref. [13| uses Uf = 5 and 



. m, and rrib, since 



A^ =321 MeV. The difference between the values of A^ is certainly insignificant. The 

MS MS ° 

different choices for m and rrib, which set the flavor thresholds, will feebly affect the evolution 
of the strong-coupling constant, the PDFs, and the FFs. In order to conservatively estimate 
the scale uncertainty, we set fiR = ^ritit, fJ'F = /u'f = ^f^t, independently vary ^r and in 
the range 1/2 < C,r,C,f < 2, and determine the maximum upward and downward deviations 
from our default predictions, for ^r = = 1. 

The theoretical predictions evaluated with the CTEQ6.5 PDFs 141 and the new FFs 

n n 

|13l | are compared with the CDF data [l| on an absolute scale in Fig. [1] and in the data- 
over-theory representation with respect to the default results in Fig. [2l The three frames in 
each figure refer to the D^, D^, and D*^ mesons; in each frame, the theoretical uncertainty 
due to the scale variation described above is indicated by the dashed lines. In all cases, we 
find that the agreement with the data is improved relative to our previous analysis jsl, as 
may be seen by comparing Figs. [1] and [2] with Figs. 1 and 2 in Ref. In fact, the 
data, except for the smallest-p^ point, and the D"^ data agree with our default predictions 
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FIG. 1: pT distributions da/dpx of p+p ^ Xc + X with (a) Xc = (b) Xc = D+ , and (c) 
Xc = D*+ for y/s = 1.96 TeV and \y\ < 1 evaluated at NLO in the GM-VFNS using the FFs 
of Ref. [13I in comparison with experimental data from CDF [1]. The solid lines represents the 
default predictions, for = S,F = ^, and the dashed lines indicate the maximum deviations for 
independent variations in the range 1/2 < S,RtCf < 2. 



within the experimental errors. This is also true for the D*~^ data, except for the two small- 
Pt points. This improvement may be attributed to the advancement in our knowledge of 
charmed- meson FFs [3] , which now also includes detailed information from the B factories, 
while the previously used FF set was solely based on LEPl data 
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FIG. 2: Same as in Figs. [Tl^a)-(c), but normalized to the default predictions. 
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FIG. 3: Same as in Fig. [2^a), but including besides the default prediction those evaluated with the 
IC parameterizations from Ref. {20I for n = 1 (solid line), 2 (dashed line), 3 (densely dotted line), 
4 (dot-dashed line), 5 (scarcely dotted line), 6 (dotted line). 
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In the following, we take the point of view that the FFs for the Xc mesons are sufficiently 
well known and ask the question whether the CDF data can discriminate between the various 
CTEQ6.5 PDF sets endowed with IC. For brevity, we only present our results for the case 
of the meson, which yields the largest cross section; the results for the and D*^ 
mesons are very similar. Specifically, we repeat the calculation of the central prediction in 
Fig. [T](a) in turn with PDF sets CTEQ6.5Cn for n = 1, . . . , 6 and normalize the outcome to 
the default prediction with zero IC of Fig. [It^a). The results are shown in Fig. [31 We observe 
that the ratios for n = 1, 2, 3, 4 lie very close to unity, the largest deviation being 2%. Only 
the ratios for n = 5, 6 significantly deviate from unity. In fact, for n = 5 the ratio is around 
1.1, and for 72 = 6 it ranges between 1.17 and 1.27 in the considered Pt range. This finding 
is easy to understand. The x values dominantly contributing to the production cross section 
at the Tevatron are typically rather small, e.g. x < xt ^ 0.025 for px = 25 GeV. In this x 
range, IC is greatly suppressed for n = 1,2,3,4, while the IC for n = 5, 6 is significant for 
X > 10^'^. The IC for n = 6, with {x)c+c = 0.024, is more pronounced than that for n = 5, 
with {x)c+c = 0.011, which explains the hierarchy of the respective ratios in Fig. [3l We 
conclude from Fig. [3] that the size of the IC enhancement is comparable to the errors on the 
CDF data only for n = 5,6. In fact, the two large-p^ data points from CDF tend to disfavor 
the sea-like IC implemented for n = 6, although a firm statement would be premature. On 
the other hand, the effects due the IC of the BHPS {n = 1,2) or meson-cloud {n = 3,4) 
models are too feeble to be resolved by the presently available CDF data. In order to obtain 
a handle on these types of IC, one needs data at considerably larger values of pr, where the 
deviations from the zero-IC prediction are sufficiently large. Of course, the cross sections 
will be much smaller there and hard to measure with sufficient accuracy. 

However, one should keep in mind that the CDF analysis of Ref. [ij] is merely based on 
5.8 pb~^ of data recorded in February and March 2002. Since then the integrated luminosity 
of run II has increased by more than a factor of 1000, exceeding 6 fb~^ as it does. Therefore, 
we conclude this section by exploring the pt range beyond 25 GeV, up to 75 GeV, assuming 
the data to be taken in the central tracking region, \y\ < 1, as before. The respective 
extensions of Figs. [T](a) and [3] are presented as Figs. Hl^a) and (b). From Fig. |4](a), we 
read off that the cross section decreases from 5 to 5 x 10"'^ nb/GeV as px runs from 25 to 
75 GeV, and that the theoretical uncertainty ranges from small to negligible. Comparing 
Fig. m^b) with Fig. [21 we observe that, in the BHPS and meson-cloud models, the sensitivity 
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FIG. 4: Same as in Figs. [IJa) and El but for 25 GeV < pr < 75 GeV. 

to IC is dramatically increased as px runs from 25 to 75 GeV, while it reaches some level of 
saturation in the sea-like model. We conclude that the measurements of the cross section 
distributions da/dpx of pp —>■ Xc + X based on the full data sample to be collected at the 
Tevatron by the end of run II would have the potential to yield useful constraints on IC. 

B. Predictions for RHIC 



Instead of measuring the cross sections at much larger values of px with the CDF and 
DO detectors at the Tevatron, it may be advantageous to lower the cm. energy to decrease 
xt in the considered pt range. A collider with smaller cm. energy is in operation, namely 
RHIC with i/i = 200 GeV, where the cross section distributions da/dpT of pp ^ X^ + X 
could be measured. So far, only the STAR Collaboration has published X^ production data, 
namely for dAn ^ + X with pt < 2.5 GeV l^]. We encourage STAR and PHENIX to 
also study pp —>■ + X at larger values of pt- 

The calculations of d^a /{dpT dy) for pp ^ Xc + X at RHIC are completely analogous to 
the Tevatron case. We only need to replace the antiproton PDFs by proton PDFs and put 
^/s = 200 GeV. We again integrate y over the range \y\ < 1. In Fig. [5l we show da/dpT as a 
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FIG. 5: Same as in Fig. [H but for pp collisions with a/s = 200 GeV. 

function of in the range 3 GeV < < 20 GeV for = D^, D+, The results for D+ 
and D*~^ are very similar, while that for has a similar shape, but differs in normalization. 
In fact, the D*^ to cross section ratio ranges from 0.46 to 0.52 in the considered px range. 
As in the case of the Tevatron, the cross sections oi pp ^ Xf. + X at RHIC rapidly fall with 
increasing value of px-, by the factor 2.4 x 10~^ as pt runs from 3 GeV to 20 GeV. Thus, 
it might be difficult to measure these cross sections in the upper part of the considered pt 
range. 

We now repeat the analysis for the meson in Fig. [5] using the various PDF sets with 
IC, CTEQ6.5Cn with n = 1,2, 3, 4, 5, 6. The results, normalized to the calculation for n = 0, 
are shown in Fig. O Here we observe a pattern familiar from the Tevatron case in Fig. El 
except that now these ratios are much larger. In fact, the results for n = 2,4 steeply rise 
with increasing value of px, reaching values of about 3 and 2.8 at px = 20 GeV. The results 
for n = 1,3 exhibit less dramatic rises because the IC is weaker in these cases. Comparing 
Figs. [3] and O we conclude that RHIC offers a much higher sensitivity to IC as implemented 
in CTEQ6.5Cn with n = 1, ... ,4 than the Tevatron. On the other hand, the results for 
n = 5, 6 at RHIC are rather similar to those at the Tevatron, being almost independent of 
Pt- But also here, the relative deviation from the zero-IC case is stronger for RHIC. 

Unfortunately, in contrast to the Tevatron, RHIC measurements at ^/s = 200 GeV will 
be limited by low luminosity. This key limitation makes the observation of high-pr Xc 
production at this RHIC energy rather unrealistic. The high-energy pp mode of RHIC, with 
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FIG. 6: Same as in Fig. [3l but for pp collisions with ^/s = 200 GeV. 

= 500 GeV, will accrue more luminosity, perhaps 200-400 pb^^. However, this will 
happen at the expense of lowering the accessible x values and, thus, of reducing the relative 
IC effects on the cross sections, especially for the BHPS and meson-cloud models. In order 
to assess the trade-off between higher luminosity and smaller relative shifts in the cross 
sections, we repeat the Z)°-meson analysis for the 500 GeV mode of RHIC. The results are 
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FIG. 7: Same as in Figs. EJa) and El but for = 500 GeV. 

shown in Figs. El^a) and (b), which should be compared with Figs. |5]^a) andO respectively. 
We observe that, as one passes from 200 GeV to 500 GeV, the cross section is increased by 
a factor of 3.6 (36) at = 3 GeV (20 GeV) in normalization, while its relative shifts due 
to IC are greatly reduced for the BHPS and meson-cloud models, as expected. 



V. CONCLUSIONS 

n 

In this paper, we updated and improved our previous analysis [5] of charmed- meson 
inclusive hadroproduction at NLO in the GM-VFNS by using as input the non-perturbative 
FFs extracted from a global analysis of e^e^ annihilation data from CESR fl])|, KEKB 12|. 



and LEPl 



24l | in the very same scheme [13|. In fact, this led to a significantly better 
description of the pt distributions of the D^, D~^, and D*^ mesons measured at the Tevatron 
as becomes evident by comparing Figs. [1] and [2] with Figs. 1 and 2 in Ref. j^. 
Encouraged by this finding, we then investigated how much room there is for the incorpo- 
ration of IC inside the colliding hadrons. Specifically, we adopted six IC parameterizations 
which are based on the BHPS meson-cloud 18], and sea-like |20j] models, imple- 



mented with two different values of (x)c+c each. For definiteness, we focused on the 
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meson. In the case of the Tevatron, we found that the BHPS and meson-cloud models yield 
insignificant deviations from the zero-IC predictions, while the shift produced by the sea-like 
model is comparable to the experimental error and tends to worsen the agreement between 
theory and experiment. However, the experimental errors in Ref. [l| are still too sizeable to 

e once the full data 

sample of run II is exploited. 

Since IC typically receives a large fraction x of momentum from the parent hadron and 
the kinematical upper bound on x scales with Xt = 2pj'/y^, the sensitivity to IC may 
be enhanced by lowering the cm. energy ^/s. This is a strong motivation for studying the 
inclusive production of charmed mesons with large values of pr in pp collisions at RHIC, 
currently operated at ^/s = 200 GeV. In fact, we found that all three IC models predict 
sizeable enhancements of the pt distributions of mesons at RHIC, by up to 75% at 
Pt = 10 GeV. We, therefore, encourage our colleagues at RHIC to perform a dedicated 
study of charmed-meson inclusive production in the high-p^ regime. In a future high-energy 
pp-collision mode of RHIC, with ^/s = 500 GeV, which is to accrue more luminosity, the 
cross sections would be increased, while the relative shifts due to IC would be considerably 
smaller for the BHPS and meson-cloud models. 
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